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Abstract. We present the three-loop QCD corrections to the quark chromomagnetic moment 
including two different nonzero masses. This is a necessary ingredient to obtain the 
corresponding corrections to the chromomagnetic coefficient in the Heavy Quark Effective 
Theory (HQET) Lagrangian. 



1. Introduction 

The anomalous magnetic moment of the electron and the muon are among the most precisely 
measured observables in particle physics. Comparing the theoretical and experimental 
predictions for the muon magnetic moment there is currently a 3.4(7 discrepancy p!] with the 
Standard Model (SM) which makes this observable very interesting at the moment. We calculate 
finite light quark mass contributions to the chromomagnetic moment of quarks, and obtain, as 
a byproduct, the corresponding corrections to the above mentioned observables and can confirm 
the results of Refs. [21 [3l [3]. Another byproduct of our calculation is the anomalous magnetic 
moment of heavy quarks, the bottom quark in particular, where we include the effect of a finite 
charm quark mass. The magnetic moment of quarks has not yet been measured experimentally, 
however, for the bottom and the lighter quarks there are upper limits from LEPl data [5]. Due 
to the lack of space in these proceedings we will present analytic results for this observable in 
Ref. [6j. 

Whereas the anomalous magnetic moments of fermions are physical observables, the 
chromomagnetic moment is not. Nevertheless, it plays a crucial role in HQET, where it 
enters the matching coefficient of the chromomagnetic interaction operator The one-loop 
correction to the chromomagnetic moment has been obtained in Refs. [3 [9]. In Refs. [101 lllj. 
the two-loop calculation has been performed, whereas light quark mass effects to this order 
have been obtained in Ref. [12]. An estimation of higher order corrections has then been given 
in Ref. [13] and the three- loop correction with one mass scale was finalized in Ref. [7]. In 
the latter the aforementioned matching coefficient of HQET is almost trivially obtained from 
the chromomagnetic moment. In the case with two mass scales, additional diagrams have to 
be calculated in the effective theory to match it to full QCD. We will present the matching 
coefficient in Ref. [6] and restrict the following discussion to the chromomagnetic moment. The 
results given in these proceedings have been published in Ref. [14j . 



Figure 1. Sample diagrams contributing to the quark chromomagnetic moment. Double 
solid lines denote light quarks, whereas solid and curly lines denote heavy quarks and gluons, 
respectively. Note that we use the background field method for the external gluon. 

2. Calculation of the chromomagnetic moment 

To calculate the chromomagnetic moment we have to consider the quark-anti-quark-gluon vertex 
in the background-field formalism in QCD. We consider the effect of a nonzero light quark mass 
at the three-loop level to this quantity. Sample diagrams which have to be calculated are depicted 
in Fig. [H When both the quark and anti-quark are on the (renormalised) mass shell and have 
physical polarisations, the vertex F^ = F^t^ can be decomposed into two form factors, 

= 7^i^i(g') - i^'y^''quF2{q^) , (1) 

where q = pi — P2 is the gluon momentum and pi and p2 are the momenta of the quark and 
anti-quark, respectively. 

The anomalous chromomagnetic moment is given by = ■^2*^-^2(0), where ^2*^ is the quark 
wave function renormalisation constant in the on-shell scheme. The total quark colour charge is 
given by Z^^Fi(O) = 1. Thus, Fi(0) is the inverse of the on-shell wave function renormalisation 
constant, which has been calculated to three-loops including light quark masses in Ref. [15j . 
Therefore, the calculation of -Fi(O) provides a strong check on the correctness of our result. 

All Feynman diagrams are generated with QGRAF |16j and the various topologies are identified 
with the help of q2e and exp [171 118j . In a next step the reduction of the various functions 
to so-called master integrals has to be achieved. For this step we use the so-called Laporta 
method |19i [20] which reduces the three-loop integrals to 27 master integrals. We use the 
implementation of Laporta's algorithm in the program Crusher [21]. It is written in C++ and 
uses GiNaC [22j for simple manipulations like taking derivatives of polynomial quantities. In the 
practical implementation of the Laporta algorithm one of the most time-consuming operations 
is the simplification of the coefficients appearing in front of the individual integrals. This task is 
performed with the help of Fermat f23] where a special interface has been used (see Ref. [24j). 
The main features of the implementation are the automated generation of the integration-by- 
parts (IBP) identities [25], a complete symmetrisation of the diagrams and the possibility to 
make use of a multiprocessor environment. As we need the form factors at zero momentum 
transfer all occurring master integrals are on-shell propagator-type integrals. They have been 
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Figure 2. Contributions of 0{e^) from the different colour structures in ([3]) as a function of 
the mass ratio x. Note the different scales of the two diagrams. 



calculated using different analytical and numerical methods, see Refs. \15\ [26] for details. To 
calculate the colour factors, we have used the program described in Ref. |27j . 

3. Results 

We write the chromomagnetic moment in the form 

n=l 



where = 0, 57721... is the Euler-Mascheroni constant, as denotes the strong coupling constant 
with rif = ni+rim + nfi active flavours (n;, rim and are the number of massless, light and heavy 
quarks, respectively). In practise, Um and nh will be equal to one, but we keep them explicitly in 
our results in order to track the various classes of diagrams in the result. We further decompose 
the three-loop contribution containing light quarks with nonzero mass into its colour structures 

jabcdjabcd^ 

+C^^^^CAnn,nm + (3) 

where Cp = {N^ — 1)/ C^^c) and Ca = -^c are the eigenvalues of the quadratic Casimir operators 
of the fundamental and adjoint representation for the SU(A'c) colour group, respectively. In the 
case of QCD we have Nc = 3 and Tp = 1/2. The dimension of the fundamental representation is 
given by Np = N,,. The symmetrised trace of four generators in the fundamental representation 
is given by df'"^df'"^ = {N^ - 1){N^ - 6N^ + 18)/(96iV2)^ present our results at the 
renormalisation scale fi = M^, where is the pole mass of the heavy quark. The nonzero pole 
parts of the various contributions in ^ are given by 
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The contribution to (7^"^*^ and (7^^*^ is presented as a series expansion up to third order in 
the quark mass ratio x = Mm/Mh, where Mm is the pole mass of the hght quark. The results 
for the finite parts of the different colour structures are given in graphical form in Fig. [2] for 
0.2 < X < 0.4, which is relevant for charm mass effects in the chromomagnetic moment of the 
bottom quark. The mass dependence of the bottom quark in the chromomagnetic moment of 
the top quark can safely be neglected and the results at x = from [7] can be used. The analytic 
results including the renormalisation scale dependence will be given in |6]. 

Acknowledgment s 

This work was supported by the DFG through SFB/TR 9 and the Graduiertenkolleg 
" Hochenergiephysik und Teilchenastrophysik" , the Science and Engineering Research Canada 
and the Alberta Ingenuity Fund. We would like to thank Jan Piclum for fruitful discussions 
about the project and Vladimir A. Smirnov for cooperation in the calculation of the diagrams. 



M. Passera, W. J. Marciano and A. Sirlin, AIP Conf. Proc. 1078 (2009) 378 arXiv: 0809. 40621 [hep-ph]]. 
R. Barbieri and E. Remiddi, Nucl. Phys. B 90 (1975) 233. 
S. Laporta, Nuovo Cim. A 106 (1993) 675. 

S. Laporta and E. Remiddi, Phys. Lett. B 301 (1993) 440.^ 

R. Escribano and E. Masse, Nucl. Phys. B 429 (1994) 19 'arXiv:he p-ph/9403304| . 
S. Bekavac, A. Grozin, D. Seidel and M. Steinhauser, in preparation 

A. G. Grozin, P. Marquard, J. H. Piclum and M. Steinhauser, Nucl. Phys. B 789 (2008) 277 ; arXiv:0707. 13881 
[hep-ph]]. 

E. Eichten and B. R. Hill, Phys. Lett. B 234 (1990) 511. 

A. F. Falk, B. Grinstein and M. E. Luke, Nucl. Phys. B 357 (1991) 185^ 

G. Amoros, M. Beneke and M. Neubert, Phys. Lett. B 401 (1997) 81 [arXiv:hep-ph/9701375|. 
A. Czarnecki and A. G. Grozin, Phys. Lett. B 405 (1997) 142 [Erratum-ibid. B 650 (2007) 447] 
arXiv:hep-ph/9701415 . 

A. I. Davydychev and A. G. Grozin, Phys. Rev. D 59 (1999) 05402 3 [arXiv:hep-ph/9809589| . 

A. G. Grozin and M. Neubert, Nucl. Phys. B 508 (1997) 311 arXiv:hep-ph/9707318] . 

S. Bekavac, PhD-Thesis, Karlsruhe, January 2009; Shaker Verlag, Aachen, May 2009. 

S. Bekavac, A. Grozin, D. Seidel and M. Steinhauser, JHEP 0710 (2007) 006 jarXiv:0708. 17291 [hep-ph]]. 

P. Nogueira, J. Comput. Phys. 105 (1993) 279. 

R. Harlander, T. Seidensticker and M. S teinhauser, Phys. Lett. B 426 (1998) 125 ;arXiv:hep-ph/9712228] . 

T. Seidensticker, arXiv:hep-ph/9905298' 

S. Laporta and E. Remiddi, Phys. Lett. B 379 (1996) 28 3 [arXiv:hep-ph/9602417| . 
S. Laporta, Int. J. Mod. Phys. A 15 (2000) 5087 |arXiv:hep-ph/0102033f 
P. Marquard and D. Seidel, unpublished. 



C. Bauer, A. Frink and R. Kreckel, [ar Xiv:cs.sc/0004015 



R. H. Lewis, Fermat's User Guide, http://www.bway.net/~lewis 

M. Tentyukov and J. A. M. Vermaseren, Comput. Phys. Commun. 176 (2007) 385 [arXiv:cs.sc/0604052| . 
K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B 192 (1981) 159. 

S. Bekavac, A. G. Grozin, D. Seidel and V. A. Smirnov, Nucl. Phys. B 819 (2009) 183 l arXiv:0903.4760l 
[hep-ph]]. 

[27] T. van Ritbergen, A. N. Schellekens and J. A. M. Vermaseren, Int. J. Mod. Phys. A 14 (1999) 41 
arXiv:hep-ph/9802376; . 



